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We consider models of heavy fermions in the strong coupUng or local moment limit and include 
phonon degrees of freedom on the conduction electrons. Due to the large mass or low coherence 
temperature of the heavy fermion state, it is shown that such a regime is dominated by vertex 
corrections which leads to the complete failure of the Migdal theorem. Even at weak electron-phonon 
couplings, binding of the conduction electrons competes with the Kondo effect and substantially 
reduces the coherence temperature, ultimately leading to the Kondo breakdown. Those results are 
obtained using a combination of the slave boson method and Migdal-Eliashberg approximation as 
well as the dynamical mean- field theory approximation. 



PACS numbers: 71.27.+a,71.10.Fd,71.38.-k,75.20.Hr 

I. INTRODUCTION AND MODELS 

The heavy fermion paramagnetic state as realized for 
instance in CeCue or induced by a weak magnetic field 
in YbRh2(Sii_a^Gea;)2 corresponds to the coherent, Bloch 
like, superposition of the individual Kondo screening 
clouds of the spins of rare-earths.^ The coherence tem- 
perature Tcoh, or inverse effective mass of this state is 
set by the Kondo^^ scale which lies orders of magnitude 
below the Fermi temperature, T/?, of the host metallic 
state. In general, for weakly correlated metals, the char- 
acteristic phonon frequencies set by the Debye tempera- 
ture Qd are much smaller than the Fermi temperature 
Ep :^ 10^~^K thus yielding a small parameter Qd/Tf 
on which Migdal theorem is based. ^ This small parame- 
ter is absent in heavy fermion materials since the Fermi 
temperature should be replaced by Tcoh which might be 
even smaller than IK. Moreover, smallness of the rele- 
vant energy scales implies that properties of the ground 
state can be easily tuned by Hamiltonian perturbations. 
This observation raises the central question of this paper: 
what role do phonons play in models of heavy fermions? 

We address this question on the basis of the periodic 
Anderson model (PAM) on a square lattice with Holstein 
phonons that couple to the conduction band electrons: 

H =Ho + i^y + Hph with 

i,(T i,a (1) 
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Here, Hq describes the conduction band with disper- 
sion relation e(fe) = — 2t(cos 

H~ cos hy) fi and with 
c\, ^ creating a conduction electron with z-component of 
spin (T and in the Bloch state with crystal momentum 
k. Next, Hy accounts for the hybridization with a con- 
duction electron in Wannier state centered around unit 
cell i and the localized /-electron in the same unit cell 
while the Coulomb repulsion set by the Hubbard U on 
the /-orbitals accounts for local moment formation. Fi- 
nally, Hph corresponds to Einstein phonons with a Hol- 
stein coupling to the conduction electrons. In fact, the 
choice between a model in which phonons couple predom- 
inately either to the conduction or /-electrons is material 
dependent. Indeed, retaining the coupling of Holstein 
phonons only to the conduction electrons is justified in 
the local moment regime where charge fluctuations on 
the /-orbitals are suppressed due to the strong Coulomb 
repulsion. In this limit the model maps onto the Kondo 
lattice model (KLM): 

Hklm = Hq -\- Hph + J ^5^5/. (2) 

i 

In contrast, the effect of phonon coupling to the /- 
electrons is expected to become more important in the 
mixed valence regime of moderately heavy fermions such 
as filled skutterudites as emphasized in Ref. [71 

We use several methods to unravel the physics con- 
tained in those model Hamiltonians. Arbitrarily low 
temperatures as well as phonon frequencies may be 
reached within a combination of the slave-boson (SB) 
mean-field approximation^ with the self-consistent calcu- 
lation of t he self -energy diagram of Fig. [T] within Migdal- 
Eliashber^^l^l^ (ME) approximation. The former leads 
to the hybridized band picture of the heavy fermion state 
while the latter allows one to account for coupling to the 
phonons. Furthermore, since the SB approximation fails 
at finite temperatures while the ME approach does not 
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FIG. 1: Self-energy diagram. The solid (wavy) line corre- 
sponds to the bare single-particle Green's function (phonon 
propagator) , respectively. 



capture vertex corrections, we equally use the dynami- 
cal mean-field theory (DMFT) approximation with a re- 
cently developed weak coupling continuous-time quan- 
tum Monte Carlo (CT-QMC) impurity solver^^ t o in- 
clude phonon degrees of freedom into the PAMfl^l^ 
Finally, in order to access the low temperature limit, 
we resort to numerical renormalization group (NRG) 
method^^ as a complementary to the CT-QMC impurity 
solver applied to the KLM 

To understand the failure of the ME approximation for 
the heavy fermion state, let us start with the phonons 
coupled to a band of conduction electrons with a flat 
density of states of width W. In this case, the self-energy 
diagram of Fig. [l] gives a mass renormalization. 



= 1 + A 1 



mo 



--^^) where A (3) 
W/2 + cjo y kW ^ ^ 



and ujQ = ^/k/M corresponds to the phonon frequency, 
and mo is the mass of the bare electron. A corresponds 
to the dimensionless electron-phonon interaction and we 
have explicitly included high frequency corrections. 

In the adiabatic limit ujq « Tcoh , we can use the above 
formula with the heavy Fermi liquid as a starting point 
rather than the bare conduction electrons to account for 
the effect of phonons. As it is reviewed in some details 
in Sec. |TI| the ratio of the coherence temperature of the 
heavy fermion state to the Fermi temperature of the host 
metal is given by the quasiparticle (QP) residue. 



I* 



(4) 



This quantity measures the overlap of the QP with a bare 
conduction electron. Hence, the coupling of the phonons 
to the heavy fermion quasiparticles will be renormalized 
by a factor Z in comparison to the coupling to the bare 
electron (i.e. g gZ). One equally expects effective 
bandwidth to be scaled by the same factor Z. Account- 
ing for these two renormalization factors, and neglecting 
the high frequency correction in Eq. ([3| gives a mass 
enhancement, 



m 



1 + AZ, 



(5) 



where m ^ ttiq/Z corresponds to the effective mass of the 
heavy fermion state in the absence of phonons. There- 



is next to unaffected by the inclusion of phonons since 
Z « 1 even in the strong electron-phonon coupling 
limit, A 1. In Sec. [ll] we verify explicitly this result 
with the use of a combination of the SB technique and 
the ME approximation. In this case the coherence tem- 
perature is also protected from adiabatic phonons which 
stems from their coupling to merely a fraction, Z, of the 
heavy QP. 

In the high frequency limit, W » uq » Tcoh the 
correct starting point is to consider the Hamiltonian 
Ho-\-Hph. Neglecting vertex corrections, one can account 
for the low energy physics of this Hamiltonian within the 
ME approximation. This approximation describes the 
formation of quasiparticles with enhanced effective mass 
set by mo(l + A) or equivalently a bandwidth reduced by 
a factor (1 + A). From there onwards one can account for 
the magnetic impurities. The hybridization matrix ele- 
ment between the QP and the /-electron is renormalized 
by a factor + A in comparison to the hybridiza- 

tion with the bare electron. Taking into account those 
two renormalization factors, bandwidth and hybridiza- 
tion, the Kondo temperature, Tk cx ex.-p{—UW/4y'^)^ re- 
mains invariant at this level of approximation. Hence, 
as in the adiabatic case, one also expects in this regime 
a very weak influence of the phonon degrees of freedom 
on the scales of the heavy fermion state. This point is 
confirmed explicitly in Sec. [IT] within the ME approxi- 
mation. 

Therefore, we anticipate that retaining only non cross- 
ing diagrams thus omitting vertex corrections leaves the 
scales of the heavy fermion state next to unaffected. Ver- 
tex corrections will clearly play a role since they lead to 
polaron binding which competes with the Kondo screen- 
ing of the impurity spins. In particular, integrating out 
the phonon degrees of freedom yields a retarded attrac- 
tive interaction. The retardation is set by I/ojq and its 
magnitude by AW/2. This term leads to binding between 
polarons into singlets and is responsible for the onset of 
superconductivity. To capture this competition and as- 
sociated substantial reduction of the coherence tempera- 
ture even at very low values of the electron phonon cou- 



pling, we carry out in Sec. HI DMFT calculations using 
two complementary impurity solvers based on: (i) a re- 
cently developed CT-QMC techniqu d^^ ' ^^ ' as well as (ii) 
zero-temperature NRG method.^ Finally, in Sec. [lV|we 
conclude and discuss the implications of our results. 

Unless stated otherwise, throughout the paper we con- 
sider the parameters U/t = 4, V/t = 1, and e/ = /i for 
the PAM and J = W/2 with W = 0.2 for the KLM, 
and choose the chemical potential /i such that the total 
particle number per unit cell reads (n) = 1.8. 
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II. 



MEAN-FIELD APPROXIMATION 



In order to understand the basic underlying physics of 
the model Eq. ([T]) it is instructive to work first in the 
mean-field framework. A good starting point to handle 
strong electron correlations is offered by the SB mean- 
field method.^ Indeed, its usefulness has been proved 
in the studies devoted to searchin g fo r the ground-state 
of both the doped Hubbard modeP^ and its multi-band 
variant.^ ^ Moreover, it has been applied to determine in- 
stabilities of the PAM,^^ and since it captures two char- 
acteristic energy scales, i.e., the coherence Tcoh and the 
Kondo Tk ones, is believed to account for the essential 
physics of the heavy-fermion systems. 

Additionally, in order to gain preliminary insight into 
the interplay between electron correlations and electron- 
phonon coupling, we combine the SB technique with 
the ME approximation.l^I^ESl Remarkably, despite obvi- 
ous flaws there are some regimes where this approach 
works quite welP^M^ and on top it might be system- 
atically improved by including the leading-order vertex 
corrections in the self-energy calculations. For example, 
it has widely been used to describe the dynamics of a 
single hole in the magnetically ordered background of 
Mott insulators and to analyze changes in the QP prop- 
erties by c ouplin g to various bosonic excitations such 
as magnons,'^^'^^ phonons,^^'^^' and orbitons.^^ Moreover, 
recent studies of the one-dimensional quarter-filled Hol- 
stein model have shown that in the temperature range 
above the onset of the Luttinger liquid phase, the ap- 
proximation is capable of reproducing the temperature 
dependence of the one-particle spectral function obtained 
within a cluster extension of the DMFT.1^ 



A. Energy scales in the PAM 
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FIG. 2: (Color online) c-electron spectral function Ac(k,c<;) 
(left) and the corresponding density of states Nc{ijo) obtained 
in the PAM with V/t = 1 in the coherent Fermi liquid regime 
at the low temperature /3t = 400. 



singly occupied /-orbitals at Tk is also seen as disconti- 
nuity of the specific heat Cy = — T( ^^^y ) shown in 
Fig. |3] The second feature of Cy is a low-temperature 
peak at Tcoh that signals entering the coherent Fermi 
liquid regime with well defined QP peaks whose effec- 
tive mass is strongly enhanced. The two energy scales 
Tcoh and Tk read off from Fig. [3] are listed in Table [T[ 
As expected based on the Gutzwiller approach,^ large 
A/'-method,EI DMFT studies,"^ as weh as on the QMC 
simulations,^ one finds that both energy scales track each 
other on varying the hybridization amplitude V. To 
be more precise within the large N approach, Tcoh lies 
well below the Kondo temperature with proportional- 
ity constant being strongly dependent on the conduction 
density of states. However, at fixed density and upon 
varying the hybridization, Tcoh tracks Tk- In particu- 
lar, for a flat density of states, the Gutzwiller method 
yields Tk cx exp(— /7W/4V^) and accordingly increasing 
V shifts both characteristic temperatures towards higher 
T. Remarkably, even though Tk changes by one order of 



It is now well established that both the PAM and its 
strong coupling version, i.e., the KLM have two param- 
agnetic solutions. The first one sets in above the scale Tk 
and is characterized by the fully localized /-electrons. In 
contrast, below the Kondo temperature Tk, the c- and 
/-electrons couple to form a heavy fermion state that 
on decreasing T evolves eventually into a coherent Fermi 
liquid state and the onset is marked by the so-called co- 
herence scale Tcoh- Its most prominent features are: (i) 
a fiat low-intensity region in the c-electron spectral func- 
tion Ac(k, cu) that results from the hybridization between 
a dispersionless correlated /-band and a wide conduc- 
tion c-band; as shown in Fig.|2]it crosses the Fermi level 
and hence it is responsible for a metallic character of the 
low-temperature ground-state; (ii) quasiparticles with a 
strongly renormalized effective mass m/mo ^ 1, and (iii) 
linear specific heat Cy = 7T. 

Within the SB approach one can define a Kondo tem- 
perature Tk- It corresponds to the temperature scale at 
which the SB factor z renormalizing the hybridization 
amplitude V vanishes. The transition into the state with 
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FIG. 3: (Color online) Temperature dependence of the specific 
heat Cv'. discontinuity of Cv in the high-temperature regime 
signals Kondo breakdown at Tk; low-temperature peaks indi- 
cate onset of the Fermi liquid regime at Tcoh- Inset shows a 
closeup of the low-temperature limit. 
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FIG. 4: (Color online) c-electron single-particle occupation 
nc(k) along the nodal direction of the Brillouin zone obtained 
in the PAM with V/t = 1. Insets show the imaginary part of 
the effective c-electron self-energy —lml2cc'^{iuJm) as a func- 
tion of Matsubara frequencies Um • 



magnitude on increasing V/t from 3/4 to 5/4, the ratio 
Tcoh/^K is next to constant (c/. Table [l|. 

The onset of the coherent Fermi liquid is also signal- 
ized by the imaginary part of the effective c-electron 
self-energy —lmT^ll^{iujm) (see Appendix). Indeed, as 
shown in the inset in Fig. [4] with V/t = 1, a conspicuous 
change in the behavior of —lmEll^{iuJrn) with ujm 
takes place in the temperature region 30 < /3t < 40 
matching the low-temperature peak in Cy. Moreover, 
a discontinuity in the c-electron single-particle occupa- 
tion nc(k) = J^^dujAc{li^uj) emerges below f3t = 100 
(see Fig. [4|. In fact, since the coherent QP part of the 
spectral function is given at low temperature by a delta 
function of weight Z: Ac(k, cj) ~ ZS{uj — £k), its magni- 
tude is given precisely by the QP residue Z which in turn 
is directly related to the effective QP mass Z~^ = m/rriQ. 
In the low temperature limit, the value of Z can be deter- 
mined as the slope of the imaginary part of the effective 
self-energy on the Matsubara axis (see the inset in Fig. [4]): 



^71 



=7rT 



Furthermore, we observe that the value of V strongly 



TABLE I: Slave-boson coherence temperature Tcoh, Kondo 
temperature Tk, their ratio, as well as QP residue Z obtained 
from Eq. ^ in the PAM at /3t = 400 for a few representative 
values of the hybridization V . 
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FIG. 5: (Color online) Low-energy part of the c-electron 
spectral function Ac(k,C(;) (left) and the corresponding den- 
sity of states Nc{u)) (right) in the PAM found at /3t = 400: 
V/t = 1.25 (top), V/t = 1 (middle), and V/t = 0.75 (bottom). 



influences the shape and the intensity of the QP band 
around the Fermi level depicted in Fig. |5j On the one 
hand, V/t = 5/4 yields a more dispersive QP band with 
a higher intensity, as compared to the V/t = 1 case. It 
extends much above the Fermi energy up to a certain 
value oJc/t — 0.06. On the other hand, a smaller V/t = 
3/4 produces a strongly renormalized low-intensity flat 
band just above the Fermi level extending up to Uc/t 
0.02. Hence, in the considered range of the hybridization 
amplitude the SB results indicate that both energy 
scales Tcoh and ujc are intimately related, i.e., Tcoh ^ 



B. Effect of elect ron-phonon coupling 

We now discuss how the previously found features of 
the coherent ground state of the PAM are altered by 
the electron-phonon coupling. At zero temperature, the 
imaginary part of the self-energy of Eq. (A. 12) reads. 



KM = -^7rwo{N+[l-e{uj + ujo)]+N-e{uj-ujo)}, 

(8) 

with = N^{uj ± uq) and 6 (a;) being a usual Heavy- 
side step function. Consequently, as shown in Fig. [6) 
^cci^) vanishes in a window —ujq < uj < ujq. Moreover, 
the hybridization gap of the heavy fermion state, equally 
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FIG. 6: (Color online) Low-energy part of the real 5]cc(^) s^^d 
imaginary — Ecc(cc;) part of the c-electron self-energy obtained 
with V/t = 1, A = 0.5, and ujo/t 0.5 at /3t 400. Inset 
depicts the corresponding c-electron density of states Nc{(jo) 
and, for comparison, Nc{uj) obtained with V/t — 3/4. 



leads to strong suppression of this quantity just above 
uo > uoq (see Fig. |6|. Away from this energy window the 
imaginary part of the self-energy is finite thereby produc- 
ing an incoherent high-energy background visible in the 
single particle spectral function (see Fig. [t]). A measure 
of this incoherent background is obtained from the value 
of nc(fc = (7r,7r)) shown in Fig. [sj Finally, since the 
imaginary part of the self energy vanishes at the Fermi 
energy, a well define QP is formed and the correspond- 
ing QP residue Z can be extracted from the real part of 
self-energy, 



J-< 



(9) 



and with the use of Eq. 

We turn now to the most important aspect of our anal- 
ysis namely the mass renormalization due to the inclusion 
of the rainbow diagrams. This quantity is plotted in Fig. 
[9]at A = 0.5 and as a function of phonon frequency uq. 
Here one finds that in comparison to the free-electron 
case, shown as inset in Fig. |9ja), the overall mass renor- 
malization is extremely small especially in the adiabatic 





(71,0) (0,0) 

k 



0.0 0.2 0.4 0.6 



FIG. 7: (Color online) Same as in Fig. [2] but with A = 0.5 
and cjo/t = 0.5. 



FIG. 8: (Color online) Same as in Fig. [2] but with A = 0.5 
and (jJo/t — 0.5. 



case. As depicted in Fig. [Toj this weak mass renormal- 
ization is also apparent in delicate flattening of the QP 
band producing very small shift in the peak position of 
the QP pole in the vicinity of = (7r,7r). 

To understand the origin of this weak mass renormal- 
ization, we can start from the hybridized band picture as 
obtained from the SB approach which up to a constant 
produces a mean-field Hamiltonian: 

HsB= ^n{k)r]l,^^r]k,a,n^Hph. (10) 

fc,cr,n=± 



where the QP energies En{k) are given by Eq. (A.7), 
Vk a n corresponding QP operators, 

while u{k) and v{k) are the coherence factors: 

1/2 



(11) 



E{k) 



1/2 



(12) 



where E{k) = y/[ef - e{k)]^ ^ 4{zV)^ 
raising and lowering operators: 



I _ cuoMQi - iPj 



Introducing the 
c^oMQi + iPi 



the phonon part of our Hamiltonian reads: 



where the phonons couple to the charge density, 



(13) 
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(15) 
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FIG. 9: (Color online) Enhancement of the effective mass 
m* obtained for a fixed electron-phonon coupling A = 0.5 
at I3t = 400 for representative values of the hybridization 
V. (a) Overall behavior of m*. Inset: enhancement of m* 
in the heavy fermion state is marginal as compared to the 
free-electron limit, (b) Low- frequency limit ujo/t < 0.2. The 
arrows indicate m* estimated from Eq. ([s}. 



In the low phonon frequency limit, it is appropriate to 
rewrite the charge density in terms of the heavy quasi- 
particles, 



{ukUk+qVk,a, + Vk+q,a,+ + ^kVk+qVk,<r,-Vk+q,a,- 
+UkVk+qVk,a,+Vk+q,a,- + ^kUk+qVk,a,-Vk+q,a, + ) ^ 



(16) 



where intra and interband transitions are apparent. For 
ujQ < Tcoh we can retain only the intra band transitions 
within the lower hybridized band, E-{k)^ thus obtaining 
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FIG. 10: (Color online) Low-energy part of the c-electron 
spectral function Ac(k,(jj) (left) and the QP peak at k = 
(tTjTt) (right) obtained with A = 0.5 and uoo/t = 0.5 at /3t = 
400: V/t = L25 (top), V/t = 1 (middle), and V/t = 0.75 (bot- 
tom). For comparison, the corresponding Ac(k = (7r,7r),c<;) 
in the absence of phonons (dashed line) is shown. 



an effective low energy Hamiltonian: 



9]J'^-^^VkVk+q4,,,_Vk+q,a,- (4 + " 

(17) 



kq,a 



This corresponds to a single band problem the band 
width being WZ where W is the bare band width, and 
with a renormalized electron-phonon coupling g gZ. 
Here, we have used the fact that in the adiabatic limit 
momentum transfer ^ ^ and Z = |v(fc)p. Hence, we 
arrive at Eq. ^ discussed in the Sec. [ij 

Let us now focus on the low-frequency limit where 
all the high-energy interband channels contributing to 
QP dressing are frozen due to the hybridization gap [see 
Fig. [9]^b)]. The sudden downturn in the data is a tem- 
perature effect. Making abstraction of this downturn by 
first taking the limit T ^ and then cjq ^ 0, gives val- 
ues of m*/m which as indicated in Fig. track the 
scale 1 + XZ thus confirming the above argument. As 
the phonon frequency grows beyond the coherence tem- 
perature, the single low energy band picture fails and in- 
terband transitions become progressively important. In- 
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deed, owing to the smallest hybridization gap (c/. inset in 
Fig. [6|, this effect appears first in the weakly hybridized 
regime with V jt = 3/4. 

Furthermore, as discussed in Sec. [T| we can acquire in- 
sight in the high frequency limit, by first taking into ac- 
count phonon degrees of freedom and then the magnetic 
impurities, to arrive at the conclusion that also in this 
limit the phonons have very little effect on the scales of 
the heavy fermion metallic state. This is again confirmed 
by Fig. [9]^a) since the overall scale of the plot remains 
very small in comparison to obtained in the absence of 
magnetic impurities. In fact, we have checked by varying 
< ojo/t < 10 and < A < 1 that the overall varia- 
tion of the Z factor remains very small. Consequently, 
we conclude that the scales of the heavy fermion metal- 
lic state are protected within the ME approximation from 
the electron-phonon interaction both in the adiabatic and 
antiadiabatic limits. This sets the stage for the impor- 
tance of vertex corrections which is the subject of Sec. 



III. DYNAMICAL MEAN FIELD 
APPROXIMATION 

We use a recently developed generalization of the weak 
coupling CT-QMC algorithm to include phonon degrees 
of freedom as impurity solver for the DMFT.^^ This 
algorithm relies on integrating out the phonon degrees of 
freedom at the expense of a retarded attractive density- 
density interaction in which one expands. This ap- 
proach to include phonons in model Hamiltonians has 
been tested extensively in the framework of the one- 
dimensional Holstein model.!^ We refer the reader to Ref. 
[T2I for a detailed discussion of the algorithm. We have 
used a stochastic analytical continuation scheme to carry 
the rotation from imaginary to real times 



pt = 200, V/t = 1, U/t = 4, <n>= 1.8, 0)=2t 



We begin with Fig. [TT] which plots the conduction elec- 
tron single-particle spectral function at our lowest tem- 
perature, pt = 200, and at the k = (tt, tt) point. As 
argued in the previous section the position of the peak 
as well as its residue is a measure of the coherence tem- 
perature Tcoh- One can estimate the QP residue directly 
from the low temperature imaginary time Green's func- 
tion by noting that 



G{k,T) 



-rA„(fc)|^^W+l|^t 



|2 for Tt » 1. 

(18) 

E^+\k) - E^+' where E^+\k) cor- 



Here, Aqp{k) 

responds to the energy eigenstate with particle number 
N -\- 1 and momentum k and Eq~^^ is the ground state 
energy in the Hilbert space with N particles. The cor- 
responding QP residue extracted by fitting the imagi- 
nary time data to the above form and in a limited range 
1 << rt < (3t/2 is shown in the inset in Fig. 
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fact that we obtain approximatively the same results for 
/3t = 120 and I3t = 200 confirms that both the tempera- 
tures lie below Tcoh of this heavy fermion state. 
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FIG. 11: (Color online) c-electron spectral function Ac{k — 
(tt, tt), uj) as a function of the electron-phonon coupling. Inset: 
QP residue as obtained from fitting the /3t = 120 (/3t = 200) 
data to the form of Eq. ( 18 ) in the range 30 < rt < 50 
(30 < rt < 70), respectively. 



As expected. Fig. [TT] tracks the evolution of Tcoh - as 
defined by the peak position in Ac {k = (tt, tt), cj) and the 
QP residue ^fc=(7r,7r) - as a function of growing electron 
phonon couplings and at fixed phonon frequency ujq = 2t. 
In the considered electron-phonon coupling range we ob- 
serve a considerable decrease of this quantity, and the 
data is consistent with a vanishing Tcoh at A 0.5. The 
fact that we cannot account for this large suppression 
of Tcoh within the ME approximation highlights the im- 
portance of vertex corrections in this parameter range. 
To confirm that the drop in Tcoh is driven by electron 



Pt=200, V/t=l, U/t=4, <n>=1.8, (0^=2t 

30 1 ^ 




FIG. 12: (a) Pair (b) conduction- and (c) /-spin suscepti- 
bilities as a function of the electron phonon coupling at our 
lowest temperature, /St — 200. 
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FIG. 13: (Color online) c-electron single-particle spectral 
function in the absence of phonons. The models parame- 
ters are given hy V/t = 1, [//t = 4, e/ = /i, and (n) = 1.8. 
The temperature from top to bottom corresponds to /St = 
10, 20, and 120. 



binding we have computed the local pairing, Xp, and spin 
susceptibilities, Xs the cluster as defined by 



X.^rd.(4(.)et(.Hc,), 



X?= [\t{S^{t)S:), 
Jo 



(19) 



where a stands for the conduction or /-electron in the 
corresponding impurity problem. Those quantities are 
plotted in Fig. [T2J As shown in Fig. 12 'a), after an ini- 
tial drop, the pair susceptibility grows and tracks the de- 
crease of the coherence temperature. The growth stems 
from conduction electron binding which originates from 
the attractive retarded interaction between conduction 
electrons mediated by a phonon exchange. Regarding a 
small drop in Xp weak electron-phonon coupling, we 



FIG. 14: (Color online) Same as in Fig. 13 but in the presence 
of phonons with A = 0.375 and ujq — 2t. 



understand this drop in terms of quasiparticle formation: 
the bare electron which enters our definition of the singlet 
5-wave pair acquires a smaller overlap with the QP as a 
function of growing A. This dressing of the bare electron 
reduces the pair susceptibility and initially overcomes the 
growth of this quantity due to electron binding. 

Electron binding into 5-wave singlets as suggested by 
'a) lead to a drop of the conduction spin suscepti- 
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Fig._ 

bility. This is clearly seen in Fig. |l2[b). Indeed, pairing 
of conduction electrons competes with the Kondo effect 
in which conduction electrons form an entangled state 
with the /-electrons thereby screening the magnetic im- 
purities. Therefore, growing values of A are linked to a 
breakdown of the Kondo singlets and hence a growth of 
the /-spin susceptibility as shown in Fig. [T2|(c). This be- 
havior should be contrasted with results obtained in com- 
plementary studies of the PAM extended by phonon cou- 
pling to the /-electrons.^ In this case, electron-phonon 
coupling generates not only pairing between /-orbitals 
but also indirectly via finite mixed valence induces bind- 
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Pt=200, V/t=l, U/t=4, <n>=1.8, X=0.25 




FIG. 15: (Color online) Single-particle spectral function of the 
conduction electrons and QP residue (inset) at fixed electron- 
phonon coupling, A = 0.25, and as a function of phonon fre- 
quency. The data points at cjq = correspond to the case 
without phonons. 



ing of the conduction electrons. 

Furthermore, Figs. [13] and [M] compare the tempera- 
ture dependence of the c-electron single-particle spectral 
functions with and without phonon degrees of freedom. 
In the absence of the electron phonon interaction, it has 
been shown^ that signatures of the hybridized bands 
in the single-particle spectral function emerge below the 
single-ion Kondo temperature. On the other hand, it is 
only below the coherence temperature that a well define 
hybridization gap is visible. As shown in Ref. [28[ this 
hybridization gap leads to a plateau feature in the mag- 
The plots in Fig. 



netization curve. The plots in Fig. |T3] do not include 
coupling to the phonon degrees of freedom and point to 
the gradual formation of the hybridization gap. It is only 
at the lowest shown temperature, = 120, that a clear 
hybridization gap is apparent. The SB mean- field anal- 
ysis gives Tx/t 0.7 for this parameter set. Hence, 
for the temperature range considered in Fig. [13] and in 
the absence of electron-phonon interaction, we always ex- 
pect signatures of the hybridized band structure when the 
data is taken below the Kondo temperature. 

Turning on the electron phonon interaction to A = 
0.375 and setting the phonon frequency to ujq = 2t re- 
duces Tcoh - as measured by the QP residue (see Fig. 
11) - by roughly an order of magnitude. Since one 
expects the Kondo temperature to track the coherence 
temperature, a similar reduction of the former quantity 
is foreseen. Our highest temperature spectral function, 
= 10, in Fig. [M] shows no sign of the hybridization 
bands thus confirming the substantial reduction of this 
scale upon inclusion of the electron-phonon coupling. At 
this temperature, only high energy features and an en- 
hancement of the effective mass by the phonon degrees 
of freedom is apparent in the data. The Fermi surface 




FIG. 16: (Color online) QP residue Z as a function of electron- 
phonon coupUng A in the KLM with J = W/2, cjq = W/2 
while W = 0.2. The conduction band filUng is Uc = 0.8. The 
red line is a guide to the eye. 



of this state accounts solely for the conduction electrons. 
Upon lowering the temperature, a hybridization gap is 
restored in the QP band but its width is substantially re- 
duced as compared to the corresponding one in Fig. [13] 
Only at our lowest considered temperature one can per- 
ceive the coherent heavy fermion metallic state with very 
faint QP peak crossing the Fermi energy in the vicinity 
of the Brillouin zone corner k — (tt, tt) thus indicating re- 
construction of a large Fermi surface including both the 
conduction and /-electrons. 

We have up to now considered only a fixed phonon fre- 
quency and altered the electron phonon coupling. Fig. 
[Ts] considers the evolution of Tcoh at a fixed electron- 
phonon interaction, A = 0.25, and as a function of fre- 
quency cjo/t. As apparent, there is a rapid convergence 
to the high-frequency limit. This crossover between the 
high-frequency and low-frequency behavior compares re- 
markably well with the energy scale of Tk as estimated 
from the SB mean- field calculations (c/. Table [l|. 

To capture the low energy behavior of the above transi- 
tion, we concentrate on the KLM of Eq. ([2]). This model 
is an effective low energy model for the PAM in the limit 
where charge fluctuations on the /-sites are negligible. 
Considering the KLM facilitates the study of the low en- 
ergy physics in this limit. We have equally used an NRG 
solveiElto at best resolve the low energies. 

In Fig. [16] we show the variation of the effective mass, 
m* oc 1/Z, induced by increasing electron-phonon cou- 
pling, A, while keeping the phonon frequency fixed at 
cjo = W 12. In analogy to our results for the PAM (see 
Fig. [TT] ) one observes a divergence of m* , or equivalently 
a vanishing of the coherence temperature. 

Next, it is particularly instructive to follow the sin- 
gle particle spectral function. Fig. [Tt] and corresponding 
density of states, Fig.[T8j as a function of electron phonon 
coupling. At weak and intermediate electron-phonon 
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FIG. 17: (Color online) Closing hybridization gap on increas- 
ing electron-phonon coupling A as seen in the c-electron single- 
particle spectral function Ac(k,cj): A = (top), A = 0.25 
(middle), and A = 1 (bottom). Parameters as in Fig. 16 




FIG. 18: (Color online) Closing hybridization gap on increas- 
ing electron-phonon coupling A as seen in the c-electron den- 
sity of states. Parameters as in Fig. [16] 



merical data which shows no sign of a coherent quasipar- 
ticle. Hence, an alternative interpretation is required. 

The characteristic feature of the state considered here 
is the frozen dynamics of the /-spin which essentially 
leads to a separation of time scales; the conduction elec- 
tron propagates in the static magnetic background of the 
/-spins. Since in the DMFT we do not allow for spin or- 
dering, it is very tempting to describe the magnetic back- 
ground in terms of an /-spin which points in a random 
direction on each site. This simple model of disorder, 
which can be dealt with within a Coherent Potential Ap- 
proximation (CPA) naturally accounts for the incoherent 
features of the low energy spectral function seen at A = 1. 
The incoherence of the single particle spectral function 
is clearly seen in the self-energy data of Fig. 



19; at A = 1 



and in the low frequency limit, the imaginary part of this 



couplings, A < 0.8, essentially the same low energy single 
particle spectral function is observed but with a renor- 
malized effective mass. Both at A = and A = 0.25 



the single particle density of states (see Fig. 18) shows 
a very clear hybridization gap which becomes somewhat 
smaller as a function of growing electron-phonon cou- 
pling thereby reflecting the reduction of the coherence 
temperature. 



In contrast, the data set at A = 1 in Figs. [T7| and 
18 stands apart. Here we are in the regime where 
the dynamics of the /-spin is frozen as signalized by 
limT^o > [see Eq. (|19|)]. In a static mean field ap- 
proach, this state is captured by the vanishing hybridiza- 
tion matrix element and thereby leads to a total decou- 
pling of the /- and c-electrons. Following this picture 
one expects the conduction-electron spectral function to 
correspond to that of tight-binding electrons coupled to 
the phonon degrees of freedom thereby producing coher- 
ent low energy single particle excitations with enhanced 
mass. However, this simple picture does not fit the nu- 
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FIG. 19: (Color online) Low frequency closeup of the imag- 
inary part of the self-energy as obtained for increasing 
electron-phonon coupling from the NRG simulations. Param- 
eters as in Fig. |16| 
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quantity remains finite!^ A similar argument was used 
in Ref.^to account for the physics of the single-particle 
spectral function in the local moment regime of the PAM 
with 5- wave superconducting BCS conduction electrons. 



IV. 



CONCLUSIONS 



In this paper, we have shown that the low coherence 
temperature characteristic of heavy fermions compounds 
leads to the breakdown of the Migdal theorem. In par- 
ticular, the dominant effect which competes with the for- 
mation of the heavy fermion metallic state is conduction 
electron binding which inhibits Kondo screening of the 
impurity spins. Since the Kondo temperature is a small 
scale, already weak electron-phonon coupling leads to a 
large reduction of the coherence temperature and ulti- 
mately to a Kondo breakdown. Our studies were con- 
ducted in the local moment regime of the periodic Ander- 
son model and Kondo lattice model in which phonon de- 
grees of freedom predominantly couple to the conduction 
electrons. Hence, we can conjecture that heavy fermion 
materials in the local moment regime should show no 
phonon anomalies since this would imply a large electron- 
phonon interaction which as we have shown leads to the 
breakdown of the heavy fermion metallic state. More- 
over, since the electron-phonon interaction considerably 
lowers the coherence and Kondo scales, it implicitly en- 
hances the importance of magnetic instabilities driven 
by the Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action between the local moments. It also leads to pairing 
hence potentially favoring superconductivity. The deli- 
cate interplay between those phonon driven instabilities 
can only be understood within the framework of cluster 
methods and is the goal of future investigations. 
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Appendix: Slave-boson technique and the 
Migdal-Eliashberg approximation 

In the SB method one enlarges the Hilbert space by 
introducing a set of four auxiliary bosonic fields e^, Pi±a, 
and di corresponding to empty, singly and doubly occu- 
pied impurity sites, respectively, such that the /-electron 



operators are replaced by /^^ fia^icn with 



= . (A.l) 



However, such enlargement of the Hilbert space intro- 
duces unphysical states which should be eliminated so as 
to recover the original Hilbert space. Therefore, the SB 
operators have to fulfill the following constraints at each 
site. 



^PlaPia = flfia- 



(A.2) 



The constraints enforced by time-dependent Lagrange 
multipliers A^^'^ and A^^'^ are conveniently handled in a 
path integral formulation. Indeed, the SB partition func- 
tion of the periodic Anderson model may be written down 
as a functional integral over coherent states of Fermi and 
Bose fields, 

ZsB = I D[e,p±^,d]Z)[c,/>[A(i),Ai')]e-(^-+'5.), 

(A.3) 

with the bosonic. 



Sb = / drJ2{4 {d^ + Af ) - Y: >^ + U)d, - Af ) 
« 

+ el(a, + AW)e, + ^pt(a,+AW-A( 

cr 

and fermionic. 



(A.4) 



P 

(A.5) 

actions. In Eq. (A.5),e/=e/ + A -^^ — /7/2 — /iisa renor- 
malized /-level energy. Next, we assume the translation 
and spin SU(2) invariance and apply the saddle-point 
approximation in which one replaces all the Bose fields 
and Lagrange multipliers by their time-independent av- 
erages, i.e., p = {p\^{t)) = (p^^(r)), and so on. The 
site-normalized SB mean-field free energy becomes. 



f3N 



^ ln(l + e-^^-^^))+/i(n) 



+ A(i)(e2 + 2p^ +d^-l)- 2A(2)(p2 + d^) + ud\ 

(A.6) 



with, 



E±{k) = - [if + e{k) ± ^[cf - e(fe)]2 + 4(zy)2j , (A.7) 

being the energies of the hybridized bands. The equi- 
librium values of the classical field amplitudes as well as 
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of the chemical potential /i are determined by the mini- 
mization FsB with respect to these parameters: 



dF. 



- / ft 



{ff)-ip' + d')=0 



(A.5 



{n)-2[{c^c) + {pf)), 



where we have defined average bond hopping Xcf — 
{c^ f)^h.c. The expectation values {a^ 13) with a = {c, /} 
are obtained from the local one-particle Green's functions 
which allows one to combine the SB approach with the 
ME approximation to account for the phonon degrees of 
freedom. The self-consistent evaluation of the diagram 
of Fig. [T] amounts to solving: 



with 



2k pN 



(A.9) 



G{k,iuJrn) 



Gcc{k, ioJrn) Gcf{k, iuJrn) 

Gfc{k,iujm) Gff{k,iujm) 

1 



(A.IO) 



Gq \k, iuj,n) 







Here D^{iVt,, 



is the bare phonon 



propagator, Um {^m) are fermionic (bosonic) Matsub- 
ara frequencies, respectively, and Go{k^iuJm) is the 2x2 
Green function matrix as obtained from the SB mean- 
field Hamiltonian at the given iteration step. Since at a 
given iteration we do not have at hand the pole struc- 
ture of Gcc{k^iuJm) in the complex frequency plane, we 
solve the above equation numerically for real frequencies. 
Namely, we use the spectral representation of the Green's 
function: 



Gcc{k,ioJm) = J duj 



;Ac{k,uj') 



where Ac{k,uj') = —^lmGlf'{k,uj'), and perform sum- 
mation over bosonic Matsubara frequencies so that the 
self-energy reads: 



^cci'^UJm) 



2k 



(A.12) 

where Nc{(jj') = -^^kAc{k^uj') is c-electron density of 
states, while n/[e(fe)] = ^^^^yj^ and 715 (cjq) = ^^n^ 
are the Fermi function and Bose-Einstein distribution, 
respectively. Hence, at a given iteration step at which 
Nc{uj) is known we compute with the above equation the 
self-energy on the real frequency axis {ioJm 0J-\-i5) and 
thereby recompute the single-particle Green's function 
and corresponding Nc{uj) until convergence of Sccl^^m) 
is reached. This in turn allows one to determine the local 
Green's functions which satisfy the usual Dyson equation 
(A.IO). Explicitly one finds: 



Gcc{kj iujm) 
Gff{k,iujm) 
Gcf{k,iujm) 



1 



■e(k) 



^eff 



^^/(k.iuj^y 



(A.13) 



zV 



^eff 



(k.iujm) ' 



where we have defined the effective self-energies: 

ZUJm - e/ 



Y^iy{k,iuJm) = [i 



- e{k) - T.^^{iuJm)\{i^ni - e/)- 

(A.14) 

Finally, the expectation values entering the SB 

saddle-point equations are obtained by performing the 
Fourier transformations: 



(at(r)/?) 



I3N 



"Gap{k,iujm), (A.15) 



and the whole process is iterated until the SB mean-field 



(A. 11) equations (A.8) are fulfilled. 
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